Search for correlations between solar flares and decay rate of radioactive nuclei 
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Abstract 

The deacay rate of three different radioactive sources ("^''K, '^^Cs and ""'Th) has been measured with Nal and Ge 
detectors. Data have been analyzed to search for possible variations in coincidence with the two strongest solar flares 
of the years 201 1 and 2012. No significant deviations from standard expectation have been observed, with a few 10"* 
sensitivity. As a consequence, we could not find any effect like that recently reported by Jenkins and Fischbach: a few 
per mil decrease in the decay rate of ^"^Mn during solar flares in December 2006. 
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1. Introduction 

In the past years, a correlation between the Sun and 
the decay rate of radioactive isotopes has been pro- 
posed. In particular, two eff'ects have been considered: 
the annual modulation due to the seasonal variation of 
the Earth-Sun distance lH] and the decrease of the de- 
cay rate during a solar flare Q. In this letter we are 
interested in the latter phenomenon. 

Briefly, solar flares are explosions on the surface of 
the Sun near sunspots. They are powered by the release 
of magnetic energy stored in the corona, up to one hun- 
dredth of the solar luminosity, and they affect all lay- 
ers of solar atmosphere, from the photosphere to the 
corona. On the Sun this amount of energy is released 
within a few minutes to tens of minutes. In this inter- 
val the plasma is heated to tens of millions of degrees 
with a strong X-ray emission and electron and proton 
acceleration (up to several tens and hundreds of MeV, 
respectively). 

In particular, the 2006 flares from December 2nd 
2006 to January 1st 2007 gave rise to X-ray fluxes 
which, measured on the the Geostationary Operational 
Environmental Satellites (GOES), were of a few times 
10"* W/m^ at the peak (see Fig. 1 of reference |2| for 
details). At that time the activity of a ~1 fiCi source of 
^*Mn was being measured by Jenkins and Fischbach Q 
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with a 2x2 inch Nal crystal detecting the 835 keV y-ray 
emitted after the electron capture decay. A significant 
dip (up to 4-10"-', ~7 o" effect), in the count rate, aver- 
aged on a time interval of 4 hours, has been observed 
in coincidence with the solar flares. On the other hand, 
a different experiment with a ~10"-^ sensitivity, carried 
out by Parkhomov |3|, did not observe any deviation in 
the activity of ''"Co, '^"Sr-Y and ^■'^Pu sources in coinci- 
dence with the same flare. 

After a few years of quiet Sun, solar activity is now 
increasing, as shown both by the increase of the steady 
X-ray flux as well as of X-flares and of other typical so- 
lar phenomena. As a matter of fact, we are approaching 
the maximum of the 1 1 year solar cycle which is pre- 
dicted to take place in Fall 2013. In our analysis we 
focus on the two most intense flares of the last years, 
namely those that occurred on August 2011 and March 
2012: X6.9 on August 9th 201 1 @ 08:08 UTC and X5.4 
on March 7th 2012 @ 00:24 UTC |4|. 

Solar flares are classified according to the power of 
the X-ray flux peak near the Earth as measured by the 
GOES- 15 geostationary satellite: X identifies the class 
of the most powerful ones, with a power at the peak 
larger than 10"* W/m^ (within the X-class there is then 
a linear scale). The two flares were well defined in time 
(a few minutes) and they illuminated the entire Earth. 
Their intensities are comparable, or even larger, than 
those observed in December 2006. During the 2011 
flare the activity of the '^^Cs and ""'Th sources were 
being measured with a Ge and with a Nal detector, re- 
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Table 1: Experimental set-ups active at the time of large solar flares. The Cs and ""'Th measurements are running underground while the K 
set-up is installed in the external laboratory at Gran Sasso. 



Flare time 


Peak intensity 


Source 


detector 


Observed decay type 


Mean counting rate 


Integration time 


August 9th 201 1 08:08 UTC 
March 7th 2012 00:24 UTC 


X6.9 
X5.4 




HPGe 


P- 


680 Hz 


1 hour 


August 9th 201 1 08:08 UTC 


X6.9 


""'Th 


Nal 


a + 13" decay chain 


3300 Hz 


24 hours 


March 7th 2012 00:24 UTC 


X5.4 


40^ 


Nal 


EC 


790 Hz 


1 hour 



spectively. On the other hand, during the 2012 flare 
the '^^Cs and "^^K sources were being studied with the 
same Ge detector and with a different Nal detector, as 
described in the next section. These different nuclides 
gives the possibility to search for possible effects corre- 
lated with solar flares in three different decay processes: 
alpha, beta and electron capture. 

2. The set-ups 

Table 1 summarizes the information on the experi- 
mental set-ups we are running to search for modulations 
in the decay rates of different radioactive sources (pe- 
riod from few days to one year). In particular, in this 
letter we only consider an interval of ±10 days around 
the time of the 2 solar flares in order to search for any 
significant deviation (positive or negative) in the decay 
rate correlated with the flares. The choice of this time 
window is quite arbitrary, since there is no model, to 
our knowledge, that coiTelates the flare intensity with 
the activity of a radioactive source. On the other hand, 
we note that, according to data shown in Fig. 2 of Q, 
the alleged influence of the flare on the source activity 
lasts for a few days around the occurrence of the flare. 

2.1. Potassium source 

A 3x3 inch Nal crystal is suiTounded by about 16 kg 
of potassium bicarbonate powder. The set-up, installed 
above ground, is shielded by at least 10 cm of lead. The 
total count rate in the 17-3400 keV energy window is 
about 800 Hz, to be compared to the background of less 
than 3 Hz when the source is removed. The energy spec- 
trum is dominated by the full energy peak at 1461 keV 
energy due to the electron-capture decay of '^''K to '**'Ar. 
The peak position and the energy resolution (^ 90 keV 
at 1461 keV) are fairly constant over months. 

2.2. Cesium source 

The activity of a 3 kBq '^'^Cs source is being mea- 
sured since June 2011. The set-up is installed in the 
low background facility STELLA (SubTErranean Low 



Level Assay) located in the underground laboratories 
of Laboratori Nazionali del Gran Sasso (LNGS). The 
detector is a p-type High Purity Germanium (96% effi- 
ciency) with the source firmly fixed to its copper end- 
cap and it is surrounded by at least 5 cm of copper 
followed by 25 cm of lead to suppress the laboratory 
gamma ray background. Finally, shielding and detector 
are housed in a polymethylmetacrilate box flushed with 
nitrogen at slight overpressure and which is working as 
an anti-radon shield. The total count rate above the 7 
keV threshold is of 680 Hz. The intrinsic background, 
i.e shielded detector without Cs source, has been mea- 
sured during a period of 70 days: thanks to the un- 
derground environment and to the detector shielding, it 
is very low, down to about 40 counts/hour above the 
threshold (0.01 Hz). The spectrum is dominated by the 
661.6 keV line due to the isomeric transition of '^^"'Ba 
from the beta decay of '^^Cs. 

Details of the experiment and the results obtained in 
the first 210 days of iTinning to search for an annual 
modulation of the '^^Cs decay constant are given in ijSj. 
Briefly, a limit of 8.5- 10"^ at 95% C.L. is set on the max- 
imum allowed amplitude independently of the phase. 

2.3. Thorium source 

The activity of a sample of natural Thorium is mea- 
sured with a 3x3 inch Nal crystal installed underground 
in the same laboratory as the Germanium experiment 
with the '^^Cs source. The sample is an optical lens, 
made by special glass heavily doped with Thorium Ox- 
ide. Note that this technique, used for improving the 
optical properties of glass, was quite common until the 
seventies. The lens is placed close to the crystal housing 
and both the lens and the Nal detector are shielded with 
at least 15 cm of lead. The total count rate above the 
threshold of 10 keV is of about 3200 Hz (gammas from 

Ac , 2i2Bi, 2i2pb, 208x1), with a background of 2.3 
Hz (due to "'°K, thorium and uranium chains and lead 
X-rays). The energy spectrum is acquired once a day, 
with a coiTected dead time of 2.63%. Even if the chain 
is not at the equilibrium, the total count rate increases 
by only 1.7 ■ lO"'* over a time period of 1 month. 
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Figure 1: Residuals, averaged over 1 day, of '"Cs and '""Th data collected around the August 9th 2011 flare and the X-ray flux measured by the 
GOES-15 Satellite. The '^^Cs data are vertically displaced by 0.001 for sake of clarity; the right-hand vertical scale gives the X-ray flux measured 
in W/m^. The two shaded bands are drawn at ±3 ■ lO"'* from the expected value. 



3. Results 

We consider separately the two largest solar flares oc- 
curred in the data taking period, i.e. X6.9 August 9th 
2011 and X5.4 March 7th 2012 [4]. For each of them 
only two of the set-ups given in Table 1 were running. 
As a matter of fact, the "'"Th set-up went out of order 
in February 2012, due to a failure in the DAQ system, 
whereas the '^''K set-up started taking data in November 
2011. On the contrary, the '^^Cs set-up is continuously 
running since June 201 1. 

Figure 1 shows the data collected in a 20 day win- 
dow centered on the August 9th 2011 flare (the day is 
given in terms of the Modified Julian Date). The X-ray 
peak flux is plotted in linear scale and given in W/m^, 
in the 0.1-0.8 nm band measured by the GOES-15 satel- 
lite m. Inside the two bands are plotted the residuals of 
the normalized count rate of the '""Th and '^^Cs sources 
(i.e. the diff'erence between the measured and expected 



count rate divided by the measured one), averaged over 
a period of 1 day. 

The error bars are purely statistical. Systematic errors 
are negligible as compared to the statistical ones during 
a data taking period of a few days only. For the ""'Th 
data a linear trend (5.7 ppm/day), due to the recover- 
ing of the secular equilibrium, is subtracted, while the 
'^^Cs data are corrected for the exponential decay of the 
source, using the nominal mean life value of 43.38 y. 
This latter correction amounts to 63 ppm/day. 

From the data we can conclude that the '^^Cs source 
does not show any significant dip or excess in corre- 
spondence with the X-ray main peak. On the other hand, 
the ""'Th source shows a questionable dip in the count 
rate, starting 1 .5 days before the X-flare. However, the 
dip is well compatible with a statistical fluctuation. As 
a matter of fact, fluctuations of the same order of mag- 
nitude can be seen at different times during the data tak- 
ing, uncorrected with X-ray flux peaks. In any case. 
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Figure 2: Residuals, averaged over 1 day, of "'Cs and ''"K data collected around the March 7th 2012 flare and the X-ray flux measured by the 
GOES-15 Satellite. The '^^Cs data are vertically displaced by 0.001 for sake of clarity; the right-hand vertical scale gives the X-ray flux measured 
in W/m^. The two shaded bands are drawn at ±3 ■ 10"'* from the expected value. 



the existence in our data of an effect as large as the one 
reported in |2|, of the order of a few per mil per day and 
lasting several days, can be excluded. The maximum 
effect compatible with our data is smaller than 3 ■ 10""* 
per day at 95% confidence level for the X6.9 flare. Such 
a limit is obtained by adding the double of the error to 
the value of the dip. 

In Figure 2 similar data for the March 7th 2012 flare 
are presented. Also in this case, no significant effect 
can be seen related to the occurrence of the X5.4 flare, 
both in the '^^Cs and in the '*"K data. An upper limit 
similar to the one given above can be issued by taking 
twice the statistical error. Note that no effect can be seen 
also in correspondence with the arrival on earth of the 
two CME (coronal mass emission) related to this earth 
facing flare, respectively on March 8th and 11th 2012 
(55994 and 55997 MJD). 

During this March 7th 2012 flare we were taking the 
'■'^Cs data also with a fully digital list mode data acqui- 



sition system, recording the time of each event. As a 
consequence, we can have the source count rate aver- 
aged over shorter time than the day or the hour. In par- 
ticular. Figure 3 shows the '^^Cs source residuals aver- 
aged over 10 minutes in a 24 hour time window contain- 
ing the X-ray peak. Again, no fast occurring effect in- 
compatible with a statistical fluctuation and larger than 
3 ■ 10"^ (at 95%C.L.) can be seen (we note that the sen- 
sitivity scales linearly with the square root of the aver- 
aging time). We cannot repeat this analysis for the 201 1 
flare because at the very timing of the flare the data ac- 
quisition has been stopped for an hour due to liquid ni- 
trogen refilling of the Germanium detector. 



4. Conclusion 

The gamma activity of three different sources, ^*'K 
(electron capture) , '^^Cs (beta decay) and ""'Th (alpha 
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Figure 3: Residuals of Cs data, averaged over 10 minutes, collected around March 7th 2012 flare and X-ray flux from the GOES-15 Satellite. 
The right-hand vertical scale gives the X-ray flux measured in W/nr. The shaded band is drawn at ±3 ■ 10"^ from the expected value. 



and beta decays) have been measured during the occur- 
rence of at least one of the two strongest solar flares 
of the years 2011 and 2012. No significant deviations 
from expectations have been observed. Up to now there 
are no quantitative models able to correlate the flare in- 
tensity with the decay constant of radioactive isotopes. 
However, from our data it is possible to conclude that 
a universal deviation of decay rate (alpha or beta or 
electron capture decay) is less than 3-10"^ per day for 
a flare of 7-10"'* W/m^ flux at the peak. By 'univer- 
sal' we mean a deviation in the count rate affecting in 
the same way all radioactive isotopes decaying through 
the same basic mechanism. We are now continuing the 
life-time measurements to search for modulations in the 
decay rate of different nuclides. This way we will have 
the opportunity to further investigate the issue of the so- 
lar flare correlation with nuclear decay in case stronger 
flares should happen closer in time to the expected 2013 
solar maximum. 
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